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ABSTRACT 

We report the discovery of a quasi-periodic oscillation (QPO) in data obtained with 
the Rossi X-ray Timing Explorer of the dipping and eclipsing low-mass X-ray binary 
EXO 0748-676. The QPO had a frequency between 0.58 and 2.44 Hz changing on 
time scales of a few days, an rms amplitude between 8% and 12%, and was detected 
in the persistent emission, during dips and during type I X-ray bursts. During one 
observation, when the count rate was a factor 2 to 3 higher than otherwise, the QPO 
was not detected. The strength of the QPO did not significantly depend on photon 
energy, and is consistent with being the same in the persistent emission, both during 
and outside the dips, and during type I X-ray bursts. Frequency shifts were observed 
during three of the four X-ray bursts. We argue that the QPO is produced by the same 
mechanism as the QPO recently found by Jonker et al. (1999) in 4U 1323-62. Although 
the exact mechanism is not clear, it is most likely related to the high inclination of both 
systems. An orbiting structure in the accretion disc that modulates the radiation from 
the central source seems the most promising mechanism. 

Subject headings: accretion, accretion disks - stars: individual (EXO 0748-676) - stars: 
neutron - X-rays: stars 
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1. Introduction 

The low-mass X-ray binary (LMXB) EXO 
0748-676 was discovered with EXOS AT by Par- 
mar et al. (1985). It showed periodic eclipses, 
irregular intensity dips, and type I X-ray bursts 
(Parmar et al. 1986). The eclipses, due to obscu- 
ration of the central X-ray source by the compan- 
ion star, occurred at a period of 3.82 hr. From 
the eclipse duration an inclination of 75° to 82° 
was derived. Dipping activity is seen at orbital 
phases ~ 0.8 - 0.2 and at (j) ^ 0.65 (0 = 
corresponds to eclipse center). The dips can be 
as deep as 100% (Church et al. 1998). It is be- 
lieved (White & Swank 1982; Walter et al. 1982) 
that such dips are caused by obscuration of the 
central X-ray source by a bulge in the outer ac- 
cretion disk that is created by the impact of the 
gas stream from the companion star, or by the 
accretion stream itself if it penetrates the disk 
further in (Frank, King, & Lasota 1987). Dur- 
ing the eclipses 4% of the 2-6 keV intensity re- 
mained, which Parmar et al. (1986) attributed to 
the presence of an accretion disk corona (ADC). 
Using their 'progressive covering' model Church 
et al. (1998) found that ~70% of the out-of-dip 
flux is contributed by this ADC, which they es- 
timated to have a radius of 0.4-1.5x10*^ cm. 

Recently, a persistent 1 Hz quasi-periodic os- 
cillation (QPO) was found with RXTE in the 
high inchnation {i ^ 80°) dipping LMXB 4U 
1323-62 (Jonker, Van der Khs, & Wijnands 1999, 
hereafter J099). J099 concluded that its pres- 
ence is most likely related to the high inclination 
of the system. In this Letter, we report the dis- 
covery of a variable 0.58-2.44 Hz QPO in EXO 
0748-676 with properties similar to those found 
in 4U 1323-62. 

2. Observation and Analysis 

We made use of data obtained with the pro- 
portional counter array (PCA; Jahoda et al. 
1996) on board the Rossi X-ray Timing Explorer 
{RXTE; Bradt, Rothschild, & Swank 1993). A 



log of all observations is given in Table |. Obser- 
vations 2 to 7 were originally performed to study 
the orbital evolution of the system. Each con- 
sisted of five ~2 ks data intervals that range in 
orbital phase from (p ~ 0.9 to ~ 0.05. All ob- 
servations yielded data in the Standard 1 and 2 
modes, which have 1/8 s time resolution in 1 en- 
ergy channel (2-60 keV) and 16 s time resolution 
in 129 channels (2-60 keV), respectively. In ad- 
dition, data were obtained in modes with at least 
2^^^ s time resolution, in 67 (obs. 1), 32 (obs. 2), 
256 (obs. 2-7), or 64 (obs. 8-10) channels, cov- 
ering the 2-60 keV range. Data during eclipses 
were removed from further analysis. X-ray bursts 
were studied separately. 

Power spectra were created in the energy bands 
2-60, 2-5, 5-8, 8-13, and 13-42 keV. The aver- 
age 1/16-128 Hz power spectra of each obser- 
vation were rms normalized, and fitted with a 
model consisting of a constant (representing the 
Poisson level), a power law, P oc z^~" (the noise 
component), and a Lorentzian (the QPO). Er- 
rors on the parameters were determined using 
Ax^ = 1. The dependence on photon energy 
of the QPO strength was determined by fitting 
the power spectra in four energy bands, while 
keeping the QPO FWHM and frequency, and the 
power law index (a) of the noise component fixed 
to the values obtained in the 2-60 keV band. Un- 
less stated otherwise the fit parameters are those 
in the 2-60 keV band. The 95% confidence up- 
per limits on the presence of the QPO and/or 
noise component were determined by fixing the 
FWHM and frequency, and/or by fixing a. 

For observations 2-6 and 8-10 we compared 
power spectra obtained in- and outside the dips. 
A count rate level was determined for each single 
orbit, below which data were assumed to be in a 
dip. For observations 8-10 this level was the per- 
sistent count rate during phases were normally 
no dipping is observed. As mentioned before, 
observations 2-6 were taken during phases that 
normally show dipping, which is confirmed by 
increases in hardness. For these observations we 
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took the post-echpse count rate as the discrimi- 
nator, since dipping usually occurs less after the 
eclipse. 

3. Results 

Figure |l] shows two typical power spectra of 
EXO 0748-676. A QPO peak was found in the 
power spectra of all observations except the first 
one. The QPO covered a frequency range of more 
than a factor 4, between 0.58±0.01 Hz (obs. 10) 
and 2.44ib0.03 Hz (obs. 2), with rms amplitudes 
between 8.4%±0.4% and 12.1%±0.6%. Within 
each observation (~1 day) the QPO was only ob- 
served in a small frequency range, never changing 
more than 15% in frequency. Table ^ gives the 
fit parameters of all observations in the 2-60 keV 
energy band. In observation 1, where the QPO 
was not detected, i.e., at least a factor 4 weaker 
than otherwise, the background corrected count 
rate was a factor 2 to 3 higher than in the other 
observations. In the other observations, neither 
the QPO nor the noise component show any cor- 
relation between their parameters and the 2-60 
keV count rate, which varies by ~20%. How- 
ever observations 8, 9, and 10, for which (due 
to the more extended phase coverage) the persis- 
tent count rate could be estimated much better 
than for the other observations, show a clear anti- 
correlation between frequency and count rate. 
The QPO peak becomes broader when its fre- 
quency increases (see Figure |[). The ratio Q, 
of frequency over FWHM of the QPO, is con- 
sistent with being constant; a linear fit gives a 
Q of 3.5±0.1, with xVd-o.f.=12.2/8. The noise 
component tends to become steeper as it gets 
stronger (see Figure The rms of the QPO de- 
creases from ~12% to ~8.5% when the rms of the 
noise component increases from ~6% to ~10.5%. 
Although there is considerable scatter around a 
linear relation, there seems to be a connection 
between the two components (correlation coeffi- 
cient is -.75). Figure ^ shows the photon energy 
spectra of the QPO in observations 2 to 10 (dip 
and non-dip data combined). It is obvious from 



this figure that they are quite flat. The energy 
spectra of the noise component are similar. 

All QPO parameters in- and outside the dips 
are consistent with each other, except for the 
QPO amplitude in observation 10, which dur- 
ing the dips is 12.0%±0.8% rms and outside the 
dips 7.5%=b0.4% rms. The noise shows only a 
significant difference in observation 4, where the 
dip rms is 9.6=b0.5% and the non-dip rms is 
14.9±1.4%. 

Four type I X-ray bursts were observed, all 
outside dips; one in observations 5 and 8 each, 
and two in observation 9 (hereafter burst 1 to 4). 
In 16 s resolution data they had peak count rates 
(2-60 keV, 3 detectors, background corrected) 
of 2330±10, 2240±10, 2410±10, and 2110±10 
counts s~^, respectively, and they all lasted -^400 
s. In all four bursts the QPO was detected. The 
QPO during burst 2 is the only that has width 
and frequency that are consistent with those of 
the QPO outside the burst. The QPO frequency 
during burst 1 is 0.1 Hz higher than that of 
the QPO outside the burst, and the FWHM 
is 0.18=b0.05 Hz, which is considerably smaller 
than 0.36 Hz. The QPOs during bursts 3 and 4 
are both poorly fitted with a Lorentzian. Both 
peaks appear to have shifted in frequency dur- 
ing the bursts. The peaked shape in the power 
spectrum of burst 4 clearly consists of two small 
peaks, which, when fitted with Lorentzians, have 
frequencies of 0.65±0.05 Hz and 0.89±0.02 Hz 
and a FWHM of ~0.1 Hz. The QPO in burst 
3 is very asymmetric. The best fit with one 
Lorentzian at 0.74±0.02 Hz and a FWHM of 
~0.1 Hz, shows excess power at higher frequen- 
cies. Dynamical power spectra of burst 3 and 4 
show that the QPO profile changed with time. 
The QPOs during the four bursts have rms am- 
plitudes of 9.7%±0.9%, 8.4%!};^%, at least 7% 
(probably 8%-10%), and ~10% (the sum of the 
rms amplitudes of the two Lorentzians), respec- 
tively, consistent with those found outside the 
bursts. 
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4. Discussion 

We have discovered a QPO in EXO 0748-676 
with a frequency varying between 0.58 and 2.44 
Hz, and with rms amplitudes between 8.4% and 
12.1%. The QPO was detected throughout the 
persistent emission, the dips, and the type I X- 
ray bursts, with strengths that were consistent 
with each other, except in observation 10, where 
the QPO was stronger during the dips than in the 
persistent emission. Over the whole frequency 
range Q was consistent with 3.5, except in the X- 
ray bursts where it was higher. The QPO was not 
detected in observation 1, when the count rate 
was a factor 2 to 3 higher than otherwise. The 
power spectral features show no changes that are 
correlated with variations of ~20% in the out-of- 
dip count rate, except for observations 8 to 10, 
in which the frequency is anti-correlated with the 
out-of-dip count rate. In addition to the QPO, a 
noise component at frequencies below 1 Hz was 
present, with a typical strength of 6%-10% rms. 
Both the QPO and the noise component have a 
flat photon energy spectrum. 

Several of the properties of the QPO, notably 
its frequency, its relatively unchanged persistence 
during bursts and dips, and its flat energy spec- 
trum, are remarkably similar to those of the ~1 
Hz QPO that was recently found by J099 in 4U 
1323-62. Figure ^ shows that the ratio of fre- 
quency to FWHM of the QPO in 4U 1323-62 
is also consistent with 3.5, and Figure ^ shows 
that the noise component in 4U 1323-62 falls 
on the strength versus slope relation found for 
the noise component in EXO 0748-676. We con- 
clude that both QPOs are produced by the same 
mechanism. Since QPOs like these have not been 
seen in other neutron star LMXBs the origin is 
most likely related to the high inclination of the 
sources (J099). It is unlikely that the QPOs are 
produced by variations in mass accretion rate 
(M) onto the neutron star: in that case they 
should also be visible in low-inclination sources. 
Moreover, other QPOs that are thought to be 



produced by M variations, show considerable en- 
ergy dependence (van der Klis 1995). A medium 
modulating the radiation from a central source 
(J099) seems to be a more promising explana- 
tion (see below). 

We have found several new features of the 
QPO, that were not seen in 4U 1323-62. First, 
the large flux range in the observations of EXO 
0748-676 allow us to say that the QPO is not 
present at high persistent count rates. In 4U 
1323-62, whose flux varied by only ~15%, the 
QPO was always present. Second, the proper- 
ties of the QPO changed during three of the 
four bursts. In 4U 1323-62 the properties of 
the QPO remained the same during all bursts. 
Third, there is the difference in frequency range 
over which the QPOs are observed: a factor 
~4 in EXO 0748-676 and ~12% in 4U 1323-62. 
This could be due to the difference in the time 
span over which the sources were observed. 4U 
1323-62 was observed for two consecutive days, 
whereas the EXO 0748-676 observations span 19 
months. The ~12% range in two days, is similar 
to the <15% range found for individual observa- 
tions (~1 day) of EXO 0748-676. Finally, the 
noise has a flat photon spectrum, which suggests 
that the mechanism behind the noise is similar 
to that for the QPO. 

Much can be learned about the QPO from 
its behavior during the X-ray bursts. The fact 
that the QPO is detected during all four bursts 
shows that the presence of the QPO does not 
directly depend on the instantaneous flux, and 
that the absence of the QPO in observation 1 is 
due to something else than just a higher count 
rate. Since the fractional rms amplitudes dur- 
ing the bursts are similar to those outside the 
bursts, variations in M onto to neutron star can 
be ruled out as the origin of the QPO. A dra- 
matic decrease in rms amplitude is expected in 
that case. The high count rates during the X-ray 
bursts enables us to detect the QPO in smaller 
time intervals, typically tens of seconds, than 
otherwise possible. In all bursts the Q of the 
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QPOs is higher than 3.5, and during three of the 
bursts we see the frequency shift. This seems to 
indicate that the QPO is intrinsically narrower 
than measured by us. The observed width of the 
QPO outside the bursts is then due to changes 
of the QPO frequency on time scales of a few 
hundred seconds, and/or due to the presence of 
several QPOs with a higher Q that occur around 
a central frequency. 

We now consider models that involve modula- 
tion of radiation from a central source. Given the 
inclination of EXO 0748-676 (75°-82°), to be in 
the line of sight to the central source, a medium 
modulating the radiation from the central source 
must at least reach a height that is between 
14% and 27% of its radial distance to the cen- 
tral source. De Jong, van Paradijs, k. Augusteijn 
(1996) found that accretion discs in LMXBs have 
opening semi-angles of ~12° as seen from the cen- 
tral source. Hence only a small height above the 
accretion disc is required, typically a few percent 
of the radial distance to the neutron star, to be 
in the line of sight. This suggests that the modu- 
lating medium is in or on the disc. Any model is 
constrained by the fact that the rms amplitude 
stays constant when the source goes into a dip; 
the medium causing the dips has to remove more 
or less the same amount of modulated as unmod- 
ulated radiation from our line of sight. Models 
can be constructed, with varying geometries that 
depend on whether the central source is point- 
like or extended, and whether the media causing 
the dips and QPO are absorbing or scattering 
media. Only a few combinations can be ruled 
out with certainty, e.g. most models involving 
a point-like source. If the dips are being caused 
by a partial covering opaque medium, the central 
source has to be extended, and the medium caus- 
ing the QPO should have an azimuthal extent 
comparable to that of the medium causing the 
dips, in order to explain the constancy of the rms 
amplitudes in- and outside the dips. The mecha- 
nism also has to produce a central frequency that 
does not vary by more than 15% over the dura- 



tion of a single observation (typically one day). 

Both partial covering by an opaque medium 
(see J099) and Thomson scattering can explain 
the flat photon spectrum of the QPO. Assuming 
the QPO frequency is the Kepler frequency of 
a structure orbiting a 1.4 Mq neutron star, one 
derives radial distances of 1.0-2.4x10® cm. The 
viscous timescale (Frank, King, & Raine 1992) 
at these radii, on which the orbiting structure is 
expected to alter, is comparable to the timescales 
that arc seen for the frequency shifts during the 
burst. The larger frequency changes between ob- 
servations might be related to the timescale on 
which the accretion disc whole changes its 
structure (typically days to weeks). We remark 
that the inferred radii are similar to those in the 
a-disc model (Frank et al.l992), at which the 
dominant opacity source changes from Kramers' 
opacity to electron scattering. This radius in- 
creases with M, and implies an anti-correlation 
between frequency and count rate, as seen in ob- 
servation 8-10. 

The fact that the QPO is not found in obser- 
vation 1 suggests that the higher count rate is 
accompanied by a change in the accretion disc 
structure. Both are probably due to an increase 
in M. Assuming that an accretion disc thickens 
with M, the modulated flux might be obscured 
by a thicker disc, but the mechanism may also 
have disappeared altogether. 

Finally we remark that the 0.15 s jitter in 
the mid-eclipse timings reported by Hertz et al. 
(1997) might be due the presence of the QPO 
reported here. A QPO with a frequency of 1 
Hz and an rms amplitude of 10%, arbitrarily su- 
perimposed on eclipse transitions lasting a few 
seconds, can be expected to cause a jitter in the 
order of 0.1 s. 
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Fig. 1. — The power spectra of observations 2 
(a) and 10 (b), showing the QPO at 2.44±0.03 
Hz and at 0.58±0.01 Hz, respectively. 
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Fig. 2.— The QPO FWHM as a function of the 
QPO frequency. The dotted hne is the best fit for 
a constant Q ((5=3.5). The black dot represents 
the value obtained for 4U 1323-62 (J099). 
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Fig. 3. — The power law index, a, of the noise 

component function of the noise rms am- 

phtude. The black dot represents the value ob- 
tained for 4U 1323-62. 



Fig. 4. — The energy dependence of the QPO 
strength (dip and non-dip data combined). The 
points without positive rms errors are upper lim- 
its. 
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Table 1 

Log of the RXTE observations of EXO 0748-676, and the power spectral fit 
parameters for each of them. the fit parameters are for all data (dip and non-dip) 

COMBINED. For EACH BURST THE QPO PARAMETERS ARE GIVEN BELOW THE OBSERVATION IN 
WHICH THEY OCCURRED; FOR BURST 4 THE FIT PARAMETERS OF BOTH LORENTZIANS ARE GIVEN. 
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d 
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PLI 
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1 


12-03-96 10 


:41 


12-03-96 


12 


;54 


6.0 (0.9) 


< 2.1= 






3.8±0.2 


0, 


,98±0, 


,12 


177.2±0.3 


2 


03-05-96 08: 


:37 


04-05-96 


00; 


;30 


12.0 (2.9) 


12 l+°'<^ 
-0.6 
8.8±0.7 


0.8±0.1 


2.44±0.03 


5.7±0.4 


0, 


,43±0, 


.09 


72.9±0.3 


3 


18-01-97 17; 


;14 


19-01-97 


13; 


;18 


9.6 (3.7) 


0.16±0.04 


0.58±0.01 


9.6±0.6 


0, 


.82±0, 


.07 


62.7±0.2 


4 


07-03-97 20; 


:03 


08-03-97 
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;40 


9.8 (5.6) 


8.9±0.5 
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0.68±0.01 


10.6±0.4 
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.06 
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6 
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q ,-1-0.8 
9.7±0.9 


0.36±0.07 
0.18±0.5 


1.30±0.02 
1.41±0.02 


7.2±0.5 


0, 


.51±0, 


.09 


64.9±0.3 


6 
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:01 
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;03 


9.7 (2.7) 
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1.04±0.02 
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0, 


.67±0, 


.08 
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7 
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;32 


8.6 


9.9±0.6 


0.31±0.05 


0.81±0.02 


S 5+" '' 

8.9±0.2 


0, 


.49±0, 


.07 


68.3±0.3 


8 


13-08-97 19 


:38 


14-08-97 


03; 


;02 


17.3 (1.7) 


9.1±0.3 


0.20±0.02 


0.833±0.005 


0, 


,55±0, 


,04 


69.0±0.1 




Burst 2 












8.4±1.1 


0.2±0.1 


0.84±0.02 
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17-08-97 12; 
Burst 3 
Burst 4 
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36.7 (6.9) 


10.4±0.2 
7.1±0.9 
7.3±1.3 
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0.11±0.07 
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0.75±0.02 
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9.7±0.2 


0, 


,54±0, 


,03 


72.5±0.1 


10 
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:13 
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18; 


:66 
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6 7+2-8 
" ' -1.6 
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0.17±0.02 


"■°°-0.05 
0.68±0.01 


10.2±0.2 


0, 


.64±0, 


.03 
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^Obs. IDs are 10068-03-01-00, 10108-01-(01-05), 20069-(02-06), and 20082-01-(01-03) for respectively obs. 1, 2, 3 to 7, and 8 to 10. 

^Time in dips (see section 2) is given between brackets. For observation 7 tliis could not be determined accurately. 

'^The average 2—60 keV count rate {3 detectors; dips, eclipses and bursts excluded), obtained from the Standard 2 data. 

'^The noise rms is integrated from 0.1 to 1.0 Hz. 

^For the frequency range 1—10 Hz and for an assumed Q of 3.5 
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